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Effect of Band Structure on the Symmetry of Superconducting States 
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Effects of the band structure on the symmetry of superconducting (SC) states are studied. For 
a square lattice system with a nearest- neighbor attractive interaction, SC states with various 
symmetries are found by changing the band structure, or, the shape of the Fermi surface. The 
spin-triplet {{px + ipy)-wave) and spin-singlet (d- or s-wave) SC states, and states with their 
coexistence {d + ipy, s + ipy) can be stabilized within the same type of interaction. The stability 
of inter layer-pairing states with line nodes is also examined, and its relation to the SC state of 
Sr2Ru04 is discussed. 
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The superconducting (SC) state of Sr2Ru04 attracts 
much attention, suicfi-|it js Ukely to have a spin-triplet 
pairing symmetryBBBQEP The triplet Cooper pairs in 
^He are formed due to the ferromagnetic spin fluctuation, 
so that it may be natural to assume that a SC state in 
Sr2Ru04 is also realized by the same mechanism. How- 
ever, recent neutron scattering experiments showed that 
the ferromagnetic {q — 0) spin fluctuation is not_large, 
but the peak is located near q = (±27r/S. ±27r/3)P con- 
sistent with band structure calculations.LP In view of this 
fact it is proposed that the antiferromagnetic spin fluctu- 
ation may lead to the spin-triplet SC statcBu-' Similar xe- 
sults have also been obtained in different contexts.ESES' 

In this article we study the effect of the band struc- 
ture on the symmetry of SC states. First we treat a 
single-band tight-binding model on a square lattice us- 
ing a mean-field approximation (MFA) . We find that the 
spin-triplet and the spin-singlet SC states, together with 
states with their coexistence can occur for the same type 
of intera&tion simply bji-changing the shape of the Fermi 
surfaceJl3) Micnas et ali3 studied the stability of various 
SC states in a model similar to ours, but they have de- 
termined only the bare T^. Namely they solved only the 
linearized self-consistency equations. We will solve the 
self-consistency equations without linearization to deter- 
mine the phase diagram, and clarify the reason for the 
change of the symmetry of the SC state as the band 
structure is changed. 

Experimental results on Sr2Ru04 seem to indicate 
that the pS.G-.^te^ in this system has a line (or lines) 
of no des .a' K4 F 9^1 [| 1- 1| ?P In order to explaiii these results 
theoreticallyji^'EBlOEj Hasegawa et al!E^ proposed an 
interlayer-pairing state which has horizontal lines of 
nodes based on a symmetry argument. We will examine 
the stability of this type of interlayer-pairing state when 
the in-plane and the interlayer interactions compete. 

First we consider a tight-binding model on a square 



lattice whose Hamiltonian is given by 
H = 



(1) 



i,3 



where ^ is the chemical potential, and is defined as 



t 



(5— ±x,ity 



s + t' 



(5— ±x±y 
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Namely, t {t') is the transfer integral for the (next) 
nearest-neighbor sites and x (y) is the unit vector in the 
X (y) direction (lattice constant is taken to be unity). 
Similarly the nearest-neighbor attractive interaction Vij 

is defined as Vij = yY.5=±x.±y^'';3+& > 0)- This 
Hamiltonian is decoupled by a standard mean-field pro- 
cedure 



(3) 



with A, 



{ciiCji) being the SC order parameter (OP). 
On the square lattice the d^2_y2- (A^), extended s- (A^), 
Px- (Ap^) andpy-wave (Ap^) symmetries are possible for 
the nearest-neighbor interaction, and the corresponding 
OP's are defined as 



Ad(z) 
A,(i) 
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Assuming that these OP's are uniform (i.e., independent 
of i) we obtain the following self-consistency equations: 



V 



i^D , > (fc)— ^ tanhf— ^ 
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where N and T are the total number of lattice sites and 



2 



Kazuhiro KuBOKl 



the temperature, respectively, and 



= — 2t(cos kx + cos ky) — At' cos k^ cos ky — ^ 



with 



j=d,S,P^,Py 
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cos fcr 



cos fcj, 
COS k„ 
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In the following we will solve the self-consistency equa- 
tions to determine the phase giagram in the plane of 
T and fi. These equations are solved by the itera- 
tion method, starting from various sets of initial values. 
When several solutions are obtained for the same set of 
parameters (/x, T), the state with the lowest free energy 
is adopted as the true one. Here we note that the SC 
long-range order cannot exist at finite temperature in a 
purely two dimensional (2D) system. However, Tc ob- 
tained within the MFA in purely 2D systems can give a 
reasonable estimate of Tc in the presence of small three 
dimensionality. 

In Fig.l we show the phase diagram in the plane of T 
and /i (or, the electron density), ioi t — 1, t' — Q and 
V = 1.5. (Due to the particle- hole symmetry, the result 
for —fj, is the same as that for ^, and so it is not shown.) 
It is seen that a (i3,2_j,2-wave SC state is stabilized near 
half- filling (/i ~ 0), while an extended s-wave state oc- 
curs at high (and low) densities {fi ~ ±4t) . In the region 
between d- and s-wave states spin-triplet (px ±zpj,)-wave 
states appear. The (px + "iPy)- and the {px — ?Pj,)-states 
are degenerate but different states, and they transform 
each other under parity V and time-reversal T transfor- 
mation. Then the system breaks V and T symmetries 
spontaneously, and these states are usually denoted as 
the chiral p-wave SC states. Near the boundary between 
triplet (pxiipy) and singlet {d or s) states we find states 
where the spin-singlet and the spin-triplet OP's coexist. 
These coexisting states, {d ± ipy)- and (s ± ipj,)-states, 
are degenerate with (d ± ipx)- and (s ± ipx)- states, re- 
spectively.cJ The (diipxHpy) and (siipxHpy) states 
are slightly higher in energy so that they are only local 
minima of the free energy. There is no reason (regard- 
ing symmetry) which precludes the coexistence of spin- 
triplet and spin-singlet SCOP's, and it is the energy that 
decides which state should appear. Actually the coexis- 
tence of d- and p-wave OP's has been found in super- 
conductor/(anti)ferromagnet bilayer systems, where the 
proximity effect induces the-imbalance of spin-up and 
spin-down electron densitiesBj-' 

The above results show that both triplet and singlet 
SC states can occur with the same type of interaction. 
In order to understand this we use the Ginzburg-Landau 
(GL) theory by expanding the free energy with respect 



to A'S 



E 



(8) 

where the gradient and higher order terms are discarded, 
and S is the area of the system. As T is decreased the 
OP with the highest (bare) transition temperature Ti"-* 
{a{Tc^^) = 0) appears. The explicit forms of a's are 
given as 



V N 



(fc) 



tanh(a/2T) 



= 2V 1 
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tanh(efe/2r) 



(9) 

The Fermi surface (FS) near the band edge {fi ^ ±4i) 

is close to the F point or fe = (±7r, ±7r), so that |ws(fc)| 

is large on the FS and thus A^ is favored. On the other 

hand the FS at half-filling is the square connecting four 

points (±7r, 0), (0, ±7r), and cosik) vanishes there. Then 

As is suppressed near half-filling and A^ is favored. For 

intermediate ji ^ ±2i), the FS comes close to the 

points kx = ±7r/2 or ky — ±7r/2 so that \u)p^^^^ \ can be 

large. Then the p-wave states have the highest in the 

region between d- and s-wave states (px and Py states 

are degenerate). 

When more than one a become negative there may 

be a coexistence of several OP's. In this case 7 terms 

will play important roles. We can explicitly show that 

P, > 0, 7|j^ > and 7^^ > {1,3 = d,s,Px.Py). The 
(2) 

fact 7,, > indicates that the OP's would form com- 
plex rather than real combinations (if they coexist), and 
in this case the nodes are removed and the system gains 
more condensation energy. This is the reason why the 
chiral [px ± ipj,)-state (rather than {px ±pj,)-statc) ap- 
pears. The coexisting states also have the complex com- 
binations of OP's due the same reason. In Fig.l the re- 
gion of ((i-|-ipj,)-state is much wider than that of (s-f-ip^)- 
state. This can be understood as follows. The nodes in 
the c?-wave state can be removed by the introduction of 
the ipx component, while in the s-state there is already 
a full gap so that the lowering of the energy due to the 
second OP is much smaller. 

Whether or not the above argument is correct can be 
tested by considering the case of t = 0, i' ^ 1, i.e., 
with only the next-nearest-neighbor hopping terms. The 
Fermi surface at half-filling consists of the lines kx — 
±7r/2 and ky = ±7r/2, so that the (p^. ±ipj,)-state should 
be most favored near half-filling. This is actually the 
case as shown in Fig. 2. Here the chiral SC state appears 
at and near half-filling, and d- or extended s-wave state 
occurs away from half-filling. 

Next we examine the stability of interlayer-pairing SC 
states in a tetragonal system with a weak interlaycr in- 
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teraction. The Hamiltonian in this case is 
H = H2D + 

i (5— ±£ i 5 

(10) 

where the summation on 5 in the second term is taken 
over 5 = ±z ± a;, ±z ± ij and iJ in eq.(l) is redefined as 
H2D- Here z denotes the unit vector in the z-direction 
with a lattice constant c. For the t_L-term we simply take 
the nearest-neighbor hopping, but we do not consider the 
nearest-neighbor interaction because of the following rea- 
son. The interaction terms with 5 = ±z (with coupling 
constant vf^^) may lead to the OP of the form A^^ij-^. 
This OP is invariant under the rotation around the z-axis 
so that it is decoupled from the in-plane OP (denoted as 
A|| ) if A|| has a d- or p-wave symmetry. Then V"^*'-' of the 
order of V (in-plane coupling constant) is necessary to 
stabilize A^^^^^. When A|| has an s-wave symmetry, it 
couples to /S.i^i±z and then the latter becomes finite even 
for an infinitesimal Since we are not interested in 

the s-wave case, we do not consider in the following. 

We decouple H^^ using the same procedure as in 2D 
case. The possible symmetries of the OP's are 

(a) sin kx cos c, sin ky cos kz c 

cos fca; sin /c^c, cosfcySin/c^c 

(7) sin/cj; sinfc^c, sin/cj^ sinfczC 

((5) cos /ca; cos fc^c, cos ky cos kzC. 

Here (a) and ((7) and (S)) are spin-triplet (spin- 
singlet) states. (Note the states with x and y inter- 
changed are degenerate.) Among these eight OP's we 
are interested in the states in (a), since their complex 
combinations, i.e., (sin kx H sin ky) coais^c are proposed 
to describe the SC state of Sr2Ru04.c3' This state has 
horizontal lines of nodes, and this behavior is consistent 
with the experimental results. For t± — all interlayer 
pairing OP's (A^) do not couple to A|| , while for t± ^ 
some of Aj_ may have a bilinear coupling (in the sense 
of GL theory) to A|| if both OP's have the same symme- 
try. In the latter case A_l can be finite once V± ^ 0. In 
view of this we take t = 1, i' = 0, Vi = 1.5 and /i = —2, 
since A|| has the {px ± zpy)-symmetry for these values 
of parameters, and OP's of (a) may couple to Aj|. In 
the following we consider (sinfc^; + i sin ky) cos k^c- and 
(sinfc^; -l- 1 sin /cj^) sin fc^c-wave OP's, and denote them as 
A||^'' and A^-*, respectively. Then we calculate An, A^-* 

(s) 

and A_^ self-consistently as functions of V±. 

In Fig. 3 the results for ij^ = at T = is shown. 
In this case A|| and A_l are not coupled. Since the FS 
has no warping along the z-axis, it is not energetically 
favorable to introduce the second component of OP (Aj^) 
because all parts of the Fermi surface are already gapped 
by A II . Thus A_L is absent for Vi_ < K (K - V). When 

V± > Vc, A^-* and A^,^'' appear simultaneously, since 
these two states are degenerate when t± . Then the state 
for large V± has (A^f' + iA^'')-symmetry and is fully 
gapped. The in-plane OP is excluded due to the same 
reason which suppresses A_l when V± < Vc- Thus we 



conclude that the state does not have nodes irrespective 
of the value of V± if t± — 0. 

Next we consider the case t± ^ 0. In this case A|| 
couples to A^''. Then the latter can be finite once V± 
becomes finite. Now the gap function is 

Afe = (sin kx + i sin fcj,)(A|| + A*]^-* cos kzC + iA^f^ sin kzc). 

(12) 

Since A*^' is induced by the bilinear coupling to A|| , their 
relative phase is either or tt, while A^'f^'' favors a phase 
±7r/2 relative to them due to 7 terms. In Fig. 4 the results 
are shown for t±_ = 0.4. Here |A|| | > |A_l| (and A^l^ = 0) 
for small V± so that there is a full gap. For large values 
of VI, lA^-*! > |A|||. However, lA^' component appears 
before |A||^^| exceeds |A|||. Then the state is again fully 
gapped except an accidental case where [A'^^^^I = \A±\. 
We have also examined other values of t±, and the SC 
state (with (px ± zpy)-symmetry for the in-plane OP) 
always has a full gap except an accidental case, unless 
t±_ becomes comparable to t. 

In summary we have studied the symmetry of the SC 
states in a single-band tight-binding model with an at- 
tractive interaction between nearest-neighbor sites. It 
is shown that the spin-triplet and the spin-singlet SC 
states, and even their coexistence can occur as the band 
structure is changed. These results can be understood by 
considering the change of the shape of the Fermi surface. 
The present result implies that the band structure is an 
important factor to determine the symmetry of the SC 
state. We have also examined the stability of interlayer- 
pairing states with line nodes. These states are difficult 
to be stabilized in a model with such a simple band struc- 
ture (Fermi surface) as that used in the present work. 
Experimental results of Sr2Ru04 seem to indicate that 
the SC state has a spin-triplet symmetry, and that there 
is a line (or lines) of nodes in the excitation gap. In order 
to give consistent interpretations of these facts, it would 
be necessary to consider the model which takes into ac- 
count thejiiore realistic crystal (and band) structure of 
Sr2Ru04.ElP This problem will be examined elsewhere. 
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Fig. 1 Phase diagram in the plane of T and /i. Param- 
eters used are t = l,t' = and V = 1.5. 

Fig. 2 Phase diagram in the plane of T and /i. Parame- 
ters used are t = 0,t' = 1 and V — 1.5. A narrow region 
between s (d) and Px + ipy states is an s -|- ipy (d + ipy) 
state. 



Fig. 3 The V± dependence of SC order parameters for 
t ^ l,t' ^ 0,V ^ 1.5, ^ = -2, T = and tj_ = 0. Note 
that all order parameters are non-dimensional. 



Fig. 4 The V± dependence of SC order parameters for 
t^l^t' ^0,V = 1.5, /.t = -2, T = and t± = 0.4. Note 
that all order parameters are non-dimensional. 
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